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Abstract: We study the exact spectrum of the AdS 4 /CFT 3 duality put forward by 
Aharony, Bergman, Jafferis and Maldacena (ABJM). We derive thermodynamic Bethe 
ansatz (TBA) equations for the planar ABJM theory, starting from "mirror" asymptotic 
Bethe equations which we conjecture. We also propose generalization of the TBA equations 
for excited states. The recently proposed Y-system is completely consistent with the TBA 
equations for a large subsector of the theory, but should be modified in general. We find 
the general asymptotic infinite length solution of the Y-system, and also several solutions 
to all wrapping orders in the strong coupling scaling limit. To make a comparison with re- 
sults obtained from string theory, we assume that the all-loop Bethe ansatz of N.G. and P. 
Vieira is the valid worldsheet theory description in the asymptotic regime. In this case we 
find complete agreement, to all orders in wrappings, between the solution of our Y-system 
and generic quasi-classical string spectrum in AdS3 x S 1 . 
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1. Introduction 

The AdS/CFT correspondence O continues to be a source of exciting new results in gauge 
and string theories. The best-studied example of the duality is the correspondence between 
four-dimensional N = 4 super Yang-Mills (SYM) theory and Type IIB superstring theory 
on AdSs x S 5 . Another example is the recently found duality between Type IIA string 
theory on AdS 4 x CP 3 and three-dimensional Af = 6 super Chern-Simons (SCS) theory p|. 
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Figure 1: Graphical representation of the Y-systems [fl9| . Circles correspond to Y- functions. Black 
nodes are the "massive" nodes which are suppressed for asymptotically large length L. For gray 
circles in the corners the equation cannot be written "locally" in terms of Y's. 



Remarkably, evidence for integrability has been found both in the gauge theory |3], |I| 
and in the string theory [||, |6| in the planar limit of large number of colors. In SYM further 
intensive development |7|, Q has led to complete description of anomalous dimensions of 
infinitely long operators by means of the Asymptotic Bethe Ansatz (ABA) equations |], [l0|] . 
Similar equations were found in [jl2| , [l3| l for SCS. Very recently the integrability approach 
was also extended to AdS3 / CFT2 dual pairs |L4|] . 

For complete solution of the planar AdS/CFT spectral problem one should be able 
to solve the integrable two dimensional worldsheet theory in finite volume. The program 
of applying the methods of relativistic integrable field theories for finite size spectrum of 
AdS/CFT was started in [15]. In [16] a generalization of the Liischer type formula was 
proposed for the first finite volume correction to the asymptotic spectrum generated by 



ABA. This information, as well as experience with relativistic integrable theories [18], led 



to the Y-system proposed in for exact solution of both AdSs x S 5 and AdS 4 x 



I>3 



theories. As we show in this work, the proposal of [19] for ABJM theory is only valid in a 



certain large subsector of the theory, and should be modified to describe the general case. 

A graphical representation of the Y-systems of |l^] is given in Fig. [I] where the Y- 
functions are represented by circles. Each value of the index A, which labels the Y-functions 
(functions of the spectral parameter u), corresponds to a node of this diagram. For each 
node A, except the gray ones, the Y-system equation has the form 

y+y- TIbV + Yb) nn 

A A ~Uc^ + yY C y (L1) 

where Y A = Ya{u ± i/2) and the index B (resp. C) labels the nodes connected to the A 
node by horizontal (resp. vertical) lines. 1 

In this paper we argue that in the AdS4/CFT3 case equation ( |1 . 1[) for black nodes 



For the gray node the equations cannot be written as functional equations in terms of y's. In many 
cases it is convenient to parameterize the Y-functions in terms of T-functions, which satisfy the Hirota 
functional equation. The "non-local" equation for the gray nodes is replaced (see for example E^|) by a 
"local" one in terms of T-functions. 
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(Fig. |], on the right) should be replaced by rather unusual equations 



(1 + 1/1^X1 + 1/1^) 

(i + 1/1^(1 + 1/1^; 



= , t ; — r — r T , a>l, (1.2) 



+ _ + _ 1 + 1® , , , 

^^-(l + i/y^)' y ^ " (1 + 1/1^) ' 



while all the other Y-system equations of [|H| need not be changed. 2 Notice that for the 
case Y >a = the new equations (|L2])-(|T~4]) coincide with the ones originally proposed in 

Once the Y-functions are found the energy of the state can be computed from 



/OO J f),mir/',.\ 
. ^7P^ ^^^~du^~^ ( f 1^^^ ) ( 1 + (^4, j ) + (^4,j ) , (1.5) 

where Uj are the exact Bethe roots given by 



YV h (u 4jj ) = -1 , Y^inj) = -1 , (1.6) 
and e is the single magnon dispersion introduced in ( [2.6D (see section 3.3 for more details). 



In the AdSs case the Y-system passes some nontrivial tests - in [19] the 4-loop pertur 



bative result [p(J was reproduced 3 , and more recently a comparison was made at 5 loops in 



22 ] . In ]23| , |24| , p6[ the Y-system was also shown to be consistent with the thermodynamic 
Bethe ansatz (TBA) approach 4 . 

The TBA equations, describing the ground state energy, do not lead to any nontrivial 
dependence of that energy on the coupling since the ground state is protected by super- 



symmetries. In [24] an extension of these equations was proposed to describe the excited 
states. These equations were solved numerically in |^9[ for the first non-trivial Konishi 
operator |30| l , giving for the first time the anomalous dimension of a non-protected operator 
in a wide range of values of the 't Hooft coupling A for a 4D gauge theory in the planar limit. 



2 With the following identification between Y-functions of [[19j and new Y-functions: Y^o = Y-* a , Ya,o = 

3 Technically the derivation of ^] is very similar to [O, where the 4-loop perturbative results were 
reproduced for the first time. 

4 In ]2(j] the Y-system was only obtained in the interval —2g<u<2g. At the same time the authors 
of failed to get the Y-system of jis| for |u| > 2g and the discrepancy was stated. The reason of this 
misunderstanding is that some of the Y-functions have branch points at u = ±2p ± | with branch cuts 
going to ±oo ± | , parallel to the real axis. Thus for real u the quantity Y(u+ |)Y(w — |) can be understood 
for instance as Y(u + \ — iO)Y(u — | + iO) or Y(u + | + iO)Y(u — ^ + iO). The first prescription (chosen in 
p^|) leads to the discrepancy whereas the second does not. Note that the problem is only present for real 
u i.e. for measure zero subset of the complex plane. Any prescription which preserves continuity leads to 
agreement with UM. In p7| (after private communication with P.Vieira) the issue was resolved. 
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The numerical results also indicate agreement with the string prediction [31] . The results 



of [^] disagree in the sub-leading 1/A 1 / 4 order with two string computations [34] and [33] 



which also disagree with each other and are based on rather strong assumptions. In [33 



a truncated model is considered whereas [34] assumes the applicability of the quasiclassics 
in the small charge limit. We hope that a first principles calculation can be done using 
Berkoviz's pure spinor formalism [37|. 



Another very recent test of the Y-system of [19| for AdSs/CFT4 was done at strong 



coupling [ 32 [ . An analytical solution of the Y-system was found for generic classical string 
motion inside AdS3 x S. It was shown to agree with the quasi-classical one- loop spectrum to 
all orders in wrapping providing thus a deep structural test of the Y-system in the regime 
where the ABA fails completely. 

In this paper we apply the technique of TBA for the AdS4 x CP 3 theory to test the 
Y-system we propose. We also present the general asymptotic infinite length solution of 
the Y-system. The asymptotic solution is very important since it allows to establish a 
correspondence between the exact solution of the Y-system and the physical states of the 
theory. It can be also used at weak coupling where it is a good approximation to study the 
leading wrapping effects. In addition, we find strong coupling solutions of the Y-system in 
two cases and compare results with the quasi-classical string spectrum thus testing deeply 
the structure of the Y-system to all orders in wrapping. 

2. Asymptotic large L solution of Y-system 

The asymptotic spectrum of the theory can be found using asymptotic Bethe ansatz (ABA) 



techniques. In this section we describe the ABA equations of [12] and link them with the 
Y-system formalism by presenting the general asymptotic solution of the Y-system. That 
solution extends the one of 

In the asymptotic regime the counting of the states is very clear and well established. 
One can analytically continue the solution of the Y-system from the asymptotic regime, 
where the solution is explicit, to finite volume. Usually this continuations is unique (see 
for example p8| ]) and allows to fix the solution of Y-system. Technically at the moment it 
is not known how to perform this procedure for the general excited state in AdS/CFT. We 



show how to apply this general method [38] for the "st(2)" subsector and also at strong 
coupling. 

2.1 Asymptotic Bethe ansatz equations for physical AdS4/CFT3 

Here we present the asymptotic Bethe equations for the AdS4/CFT3 theory, which were for 



the first time obtained in [12]. We will also introduce some notation useful for the sequel. 



J Much later the equations of [|29|] were rederived by another group |58[. The authors of |5q] confirmed the 
validity of the equations at least in the range of the coupling < A < 700 where a numerical solution was 
obtained. The perturbation theory for the world-sheet sigma model in the formulations of pvj] is naturally 
organized in powers of 1/vA and thus the value A ~ 700 should already give the asymptotic of E(X) with 
a good precision especially when an appropriate extrapolation procedure is applied. This holds assuming 
the analyticity of E(\) for real positive values of A which is however doubted in ]5q]. 
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First we define the Zhukowski variable x(u): 



1 u , 

where h(X) is some unknown function of the 't Hooft coupling A. It should have the 
following asymptotics at weak coupling and strong coupling: 

h{\) = A + h 3 X 3 + O (A 5 ) = ^/\/2 + h° + O ( -J= ) • (2.2) 



Vx 

Recently the coefficient /13 = —8 + 2£j was computed directly from the Super-Chern- 
Simons perturbation theory ||. At strong coupling the situation is less clear: in [^] and 
1 39] the coefficient h° was argued to be whereas in |H| some evidence was given in favor 



of a different value — (see also [35]). Hopefully this issue could be analyzed from 
world-sheet sigma model first principles calculation like in [4C]. 



Equation ([2JJ) admits two solutions, and we define two branches of the function x(u), 
which are called "mirror" and "physical": 



Here, by y/u we denote the principal branch of the square root. This definition of mirror 



and physical branches is the same as in the AdSs / CFT4 case |2T], |2J] , with the AdSs / CFT4 
coupling g replaced by h{X). Above the real axis, the mirror and physical branches coincide. 
x ph (u) is obtained by analytical continuation from the upper half plane to the plane with 
the cut (— 2h, 2h), and x miv (u) - by continuation to the plane with the cut (— 00,— 2h) U 
(2h, +00). The AdS4/CFT3 Bethe equations [12] for the original (physical) theory are 



written in terms of x ph (u), while the mirror Bethe equations we conjecture include x mn (u), 
in analogy with the AdSs/CFT4 case |2l| (see section 3). In sections 4 and 5 we use the 
mirror branch of x if its argument is a free variable, and the physical branch for x(uj), 
with Uj being the Bethe roots. 

In the physical ABA equations of []12| there are five types of Bethe roots: ui, U2, U3, U4 
and ug. Conserved local charges (the heights Hamiltonians) in AdS4/CFT3 are expressed 
in terms of the momentum-carrying roots 114 and u^. 

Qn = Yl <lr>4,,) + Yl <ln(ui,j) , q„ = -^—r ( )n _ x ~ ) , (2.4) 

j=l j=l 

where we have used general notation 

f ± (u) = f(u±i/2), /M = /(« + io/2). (2.5) 

In particular, string state energies in AdS4 x CP 3 or operator anomalous dimensions in the 
dual gauge theory are obtained from E = h(X)Q2- 
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The momentum and energy which correspond to a single Bethe root 114 or ug are given 



by 



% - x- 2 
and the charge Qi is the sum of all momenta: 

K 4 K. 



i=i i=i 



(2.6) 



(2.7) 



To write the Bethe equations in compact form, we introduce the following notation: 

Ki 

ti . i \ -i.j D _ t 



11 CtTU/2 
1=1 



J 



Rl = \\ix{u) - Xlj) 
i=i 



H (±) = TT l l x ( u ) x l,j 
j=l t X l,j> 



K, 



Ki 

j'=l 

A', 



3, 



(2.8) 

(2.9) 
(2.10) 



i=i 



i=i 



where cbes is the Beisert-Eden-Staudacher dressing kernel [12|. The Bethe equations of [12 
in 5(2 favored grading have the form 6 



+1 



n+nHoH 
-\iQi ^2-^4 ^4 



Q2 Q1Q3 



«i,fe 



«2,fc 



+1 = e5 



Q2 Qi Q3 



i iQl Qt R i ]h } 



Q2 R4 -R| 



(2.11) 



«3,fc 



-+ D -(-) D -H 



■ 1 _ P |<ei e -wp(«4,*) B i R 3Q4 Rj R4 

+ i ~ 6 6 +(+) +(+)° 4 ° 4 



-SaSa 



B 1 R 3 Q i R 4 Rj 



+1 



e 2«Sl e -ijp(M4 jfe ) 



B\ R3 Qi R±^R^ + 



a 4.k 



«4.fc 



where L is the length of the effective spin chain, and corresponds to the string momentum 
or length of the operator in the CS theory. The above equations describe the spectrum 
correctly in the limit L — > 00. We stress again that in those equations the physical branch 
x ph of the function x should be used in all places, e.g. inside expressions ( fDj|) , ( [2.9D , (|2.1C| ) 
for Bi, Ri and 5*/. 



6 Here, as well as when constructing the asymptotic solution of Y-system, one should be careful with the 
sign ambiguity in the square root factors inside e^* 21 and Bi,Ri. 
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The Bethe roots are additionally constrained by the zero momentum condition 



K 4 T + K- 4 + 

i= n— n^^Qi= 2 ™- ( 2 - 12 ) 

i=i X4 J i=i x i,j 
2.2 General asymptotic solution 

As we mentioned in the beginning of this section the asymptotic (large L) solution of the 
Y-system plays an important role in the whole Y-system construction. It allows to link 
a particular solution of the Y-system with an actual state of the theory. The asymptotic 
solutions are in one-to-one correspondence with the solutions of ABA equations. 

In many cases one can analytically continue a solution from asymptotically large vol- 
ume to finite volume. In J3^] another way to inject information about the state of the 
theory was proposed: demanding that the exact functions Y as approach the formal asymp- 
totic solution for infinite a or s 7 . Then one can still use the same counting of the states 
as in the ABA even for finite volumes 8 . This prescription was shown to work especially 
successfully in the strong coupling scaling limit ||32||, which we describe below. 



In view of its importance we will review the construction of [19] for the asymptotic 
large L solution of the Y-system in this section and extend it to the case ^ Y+. To 
distinguish the asymptotic Y functions from the exact ones we use the bold font: 



rp+ rp— T+ T 

- L a,l- L a,l -, - L l,s- L l,s 



1 , l/Y Gs = - ^ ^ - 1 (2.13) 

J-a+1,1 J-a-1,1 1 l,s+l 1 l,s-l 
/ \ L a ~ 1 

Y ^^ ™ T «M II n na tu + in)$l™(u + in) , (2.14) 




a-1 
2 



a-1 
2 



Y *"-™ T °' 1 II ^r^ + ^Wr (u + in) (2.15) 



a-1 
2 



where 6^ a is for even and 1 for odd terms in the product: 



E _ J 1, n + is even 
na = I 0, n + s=i is odd 



(2.16) 



and = 1 — @na- The factors &i(u) and $g(u) are constructed in such a way that the 
ABA equations fl2.ll ) for the momentum carrying nodes are given by Y£(u4j) = — 1 and 
Y(P h (u2j) = -1. This leads to (using that Ti^uzy) = -Qt/Qz) 

Ba} ^ R-2 ^ B^ Br> m B^~^ R^i ^ B^ Bo , /r, 

*4(«) = Stfi * * 1 * e-*W , * 4 (u) = S 4 S- 4 y- * 1 3 e+* gl/2 ■ 



B^-R^B^Bt ' 4W 4 BW"4 +)+ Srfl3 



(2.17) 



7 This should give the same result as analytical continuation in L. Usually, variations of Y's in L vanish 
at large a and s. The values of the Bethe roots inside the asymptotic solution should be equal to their 
exact values e.g. I^ h (u4,j) = — 1. One should study this point in more detail. 

8 One cannot exclude completely that this procedure fails for some particular small volumes. 
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The Ta.s functions which enter the definitions of Y QiS can be computed from the generating 
functional g|, @3] 



q^+)+rW- 

Q+B(-)+R(-)~ 



Q+QJ-fl(+)- ^ 



^ 3 V2 



R(-)- 



where -D 



(2.18) 

is the shift operator and R = R4R4, B = B4B4. Expansion of this 
generating functional yields eigenvalues of the su(2|2) transfer matrices: 

00 00 
W = Y,TiAu + i^)D s , W- 1 = Y,{-V a TaAu + i^)D a . (2.19) 



s=0 



a=Q 



In Appendix B we also present the expressions for the asymptotic solution after the 
duality transformation, which exchanges the sl(2) and su(2) sectors. One can see from 
those formulas that for u^j = ug j the asymptotic solution exactly conicides with the one 
proposed in IjUfl . 

In the next subsection we expand the asymptotic solution in the scaling strong coupling 
limit. 

2.3 Asymptotic solution in scaling limit 

The scaling limit is the strong coupling limit A — > 00 where the number of Bethe roots 
M and the operator length L go to infinity as vA. The Bethe roots Xj are distributed 
along cuts C on the complex plane x in this limit [^TJ. These cuts can be understood as 
branch cuts of a 10-sheet Riemann surface which corresponds to a certain function. One 
can interpret them as the eigenvalues of the classical monodromy matrix, which are usually 



written as A a = e iqa , with q a being the so called quasi- momenta. Similarly to [ 12 1 for the 
77 = — 1 grading we get 

Lx/h+Q2X 



Q2 

q 4 

Qa 



Lx/h- 



Lx/h- 



-1 



+H\ — H4 — Hg + H3 
— H2 + H\ + — H2 
—Hz + H2 + H2 — Hi 
H3 — Hi 
+H4 — H^ + H4 — Hi 
a = 6 10 , 



(2.20) 



where the resolvents H a have the form 

,.2 



H a (x) 



1 



1 X %a,j 



H a (x) = H a (l/x) . 



In these terms the Bethe equations ( [2.11 ) are equivalent to the condition that the two 
eigenvalues of the monodromy matrix are equal along the branch cut 



qi{x + iO) — qj(x — iO) = 2nn , x € C . 



(2.21) 
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We can now simplify (|2.18| ) for strong coupling. First of all we notice that the shift 
operator D becomes a formal expansion parameter. Then we use 



q+bH+rFP 

QjQ 2 -R^- 
e¥^Q~Q 2 R{-)- 

el^Q-Q 2 i?(-)- 
Qt 



exp [-i {H x - H 4 - H 4 + H 1 + H 4 + H- 4 )] 
Qi+xQ 2 



~ cxp 



~ exp 



x 2 — 1 

gi + xQ 2 
x 2 — 1 



+ Ho — Hs — H 3 + i?o ~~ -f^4 — -f^2 



+ -f^2 ~~ Hi — Hi + — H 4 — Hi 



~ exp[+i(# 3 + i? 3 )] . 



The generating functional ( [2.18 ) becomes 

w= (1-Aid)(l-A 2 d) 
(1-A 3 d)(l-A 4 d) ' 

where we have redefined the formal expansion parameter in the following way 

/ Lx/h + xQ 2 



d = exp 



+ H 4 + H 4 -H 1 + H, 



D 



Expanding the generating function ( 2.23 ) we get 

A 4 S - 1 (A 4 - Ai)(A 4 - A 2 ) - A 3 S - 1 (A 3 - Ai)(A 3 - A 2 ) / d 



Ti,s 



A 4 — A 3 



D 



a Xi a ~HXi - A 3 )(Ai - A 4 ) - A 2 a - 1 (A 2 - A 3 )(A 2 - A 4 ) / d \ a 



Ai — A5 



D 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



It is now straightforward to compute Y& a and Yq, from ( |2.13| ). Note that the factors (-^)' 
and (-^) a are irrelevant here and thus Ya are rational functions of A a only! 
Moreover, using the relation 



— - ) $4(u) ~ exp 



/{ xL/h + X Q 2 +2H __ S4 + B B B 
x z — 1 



(2.26) 



and the same relation with 4 and 4 exchanged, we obtain expressions for the massive nodes: 



v 1 \ ^1" x (Ai — A 3 )(Ai — A 4 ) — A 2 a 1 (\ 2 — A 3 )(A 2 — A 4 ) 
Y -a - (-*5) 7 7 , 

Ai — A 2 

v ( \ \+"a Ai a_1 (Ai - A 3 )(Ai - A 4 ) - A 2 a_1 (A 2 - A 3 )(A 2 - A 4 ) 
Y *-a — v~ A 5j 7 7 , 

Ai — A 2 

which are again written solely in terms of the eigenvalues of the classical monodromy 
matrix! Here, we have introduced u a , which is defined to be 1 for odd a and zero for even 
a. 
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3. TBA equations for AdS 4 /CFT 3 



In this section we derive the Thermodynamic Bethe ansatz equations for AdS4/CFT3. Let 
us first describe the general form of the TBA method [57] (see a nice introductory paper 



J21[|). We start with an integrable quantum field theory in 1+1 dimensions, on a circle of 
circumference L. The partition function of this theory at temperature 1/R is 

Z(L,R) = Y^e~ REkiL \ (3.1) 

k 

and in the limit R — > oo, R ^> L we have 

Z{L,R) ~ e - REo(L \ (3.2) 

where Eq(L) is the ground state energy. Denoting by <p and ip the bosonic and fermionic 
fields, respectively, we can write the partition function as a functional integral 

Z(L, R) = J V<pV^e' SE (3.3) 

where Se is the theory's Euclidean action. In this integral, fermionic fields are periodic 
(resp. antiperiodic) in space (resp. time), while bosonic fields are periodic in both space 
and time: 

ip(x + L,t) = tp(x,t), ip(x,t + R) = -ip(x,t) 
<j)(x + L,t) = <t>(x,t), 4>(x,t + R) = 4>(x,t). 

Using this representation of the partition function, one can relate it to the Witten index 
of the "mirror" theory in volume R: 

W(R, L) = Y^(-l) F e- LE " I{R) = £ e- RE *M = Z (L, R) . (3.5) 

k k 

The mirror theory is obtained from the original one by a double Wick rotation, and F 
in (|3.5| ) is 1 for fermionic states and otherwise. Introducing the mirror bulk free energy 
J rmir (L), defined by the mirror theory's Witten index at temperature 1/L, 

-RLF miT (L) = ]nW(R,L), (3.6) 

we see that the finite volume ground state energy is related to the infinite volume mirror 
free energy: 

Eo(L) = LF miv (L). (3.7) 

The mirror theory's infinite volume spectrum is described by the ABA equations, which 
allow one to find J 7 ™ 11 ^) and then the original theory's ground state energy. 

To compute J 7 ™ 11, it is essential to know the structure of the solutions of infinite volume 
mirror ABA equations. For numerous theories (see [|l8|), so-called string hypotheses have 
been formulated, which describe the complexes Bethe roots form in the infinite volume 
limit (simplest of those complexes are strings of roots). We will use indices A, B, ... to 
label the complexes, and denote the energy and momentum of a complex by, respectively, 
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ip* A and ie A , to underline that the mirror theory is obtained from the physical one by a 
double Wick rotation. 

Multiplying the Bethe equations for all roots in a complex, one obtains equations 
for the density pa(u) of complexes, with u 6 M. being the center of the complex. Those 
equations have the form 

Pa{u) + PA {u) = _ K BA (v, u) * PB (v) (3.8) 

where p is the density of holes, K{v,u) * f(v) = J_°° dvK(v,u)f{v) and summation over 
B is assumed. Also, we use the normalization 



/ 



, . , total number of complexes of type A 

dup A {u) = . (3.9) 

K 



The free energy is given by the minimal value of a functional of the densities 

F™{L) = min du (i Li P*A + h A )p A - PA log (l + ^ J + p A log (l + ^ J V 

(3.10) 

with constraints ( |3.8|) on the densities. Here, Ha = log[(— 1) Na ] , where Na is the number 
of fermionic Bethe roots in the complex A. Minimization of this functional gives the TBA 
equations 

lo g y A (u) = K AB (u,v) * log[l + l/y B (v)] + iLp* A + h A , (3.11) 

where y a = and K(u,v) * f(v) = f dvK(u,v)f(v). Lastly, the free energy can be 
expressed in terms of a solution of TBA equations: 

A J 

From the free energy, one can compute the ground state energy of the physical theory 
via (|3.7|). In addition, the TBA equations can be modified in such a way that their solutions 
provide also energies of certain excited states in finite volume. 

3.1 Ground state TBA equations for AdS 4 /CFT 3 . 

We first present, as a conjecture, the ABA equations for the mirror of AdS4/CFT3 theory. 
Like the physical Bethe equations [12], those equations involve Bethe roots m, «2, «3, «4, «j, 
with all roots except U2 being fermionic. The only roots which carry energy or momentum 
are 114 and u%. For a single root, we denote energy by ip\, and momentum by ie J , where 

Pi = -log—, 4 = l + h(X)(^-^). (3.13) 
ix 2, \x + x J 

Here, and everywhere in section 3 unless otherwise stated, we use the mirror branch of the 
function x(u). Note that pi (resp. e*), evaluated in physical instead of mirror kinematics, 
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coincides with the momentum (resp. energy) of a single Bethe root in the physical the- 
ory |ji~2| . This is in accordance with the fact that the mirror theory is obtained from the 
physical one by a double Wick rotation. The momentum/energy in mirror and physical 
AdSs/CFT4 are related in a similar way. 

The mirror Bethe equations we propose are written in terms of the functions B[, Ri, Si,Qi, 
which were introduced in section 2 (note that x in them should now be understood as x mir ). 
The equations for m,v,2 and 1*3 are: 



1 



Q 2 B 4 B- 



(+) D (+) 



1 



Q2 QiQ?j 



«i.fe 



q 2 r 4 r\ ^ 



Q 2 R^ R 



(+) „(+) 



(3.14) 



u 3,k 



The r.h.s. of the equations for u\ and U3 is not always unimodular, because £?, («) and 
R^\u) have cuts on the real axis. However, in the thermodynamic limit (see below) 
the single fermion roots U\,Us are distributed [26] within the interval — 2h < u\ < 2h, 
—2h < 113 < 2h, and unimodularity of the r.h.s then follows. Note that no conditions have 
to be imposed on the 141,143 roots which are parts of pyramid complexes A n , as the terms 
containing cuts cancel during fusion of Bethe equations. This can be seen from the fact 
that the kernels JC(u,v) in TBA equations (see below) for interactions involving pyramids 
are real for real u, v and have no cuts on the real axis. 
The equations for momentum-carrying roots are: 



3 i?e*(«4 jfc ) 



BfR^Q 4 + R 4 ^ ^R- 4 I - 1 
B{ R' 3 Q 4 ^R 4 ^R^ + \ 



S 4 S 4 



4,fc . 



K 4 



n 



K- A 



11 



1 \ X 4j j=1 \ 



2d. 
(3.15) 



for u = u^k-, and for u = we have 



(-) 



B+R+Qt + R, ( - ] R 



BiRlQ, 



(+) 



Note that the combination 



a 



3=1 



{x(u),xi tj ) = n \ 
3=1 \ 



S 4 S 4 



x i,j x ~ 

X T,3 X ~ 



4,fc 



n 



1 \ X 4,j j=1 \ 



X 



X 



M 

2,i 
(3.16) 



Ki 



■<TBEs{x(u),Xij) = Sl(u)Y[ 



X tf_ 



(3.17) 



is a unimodular function (see [24|, |25[| ). By ctbes(x(u),xij) we denote the usual Beisert- 
Eden-Staudacher dressing kernel analytically continued from Imit > i/2 between the 
branch points u = ±2h + i/2. 9 The above equations are similar to the Bethe equations 
for physical ABA ( |2.11| ). The difference is in the choice of the mirror branch of x(u), 
interchange of the energy and momentum (with multiplication by i) and various factors of 
\J x + jx~ , tuned in such a way that the right-hand sides are unimodular functions. This 
prescription is based on the corresponding conjecture in the AdSs/CFT4 case [ElL 



9 "Physical" choice of the branch corresponds to analytical continuation to the plane with the cut 
i/2, 2h + i/2]. In the "mirror" kinematics all cuts should go through infinity. 



-2h + 
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• 


U4 


• 





Odd-Even complexes ■% and ^4 Even-Odd complexes and +4 

Figure 2: Strings of alternating roots on the {ReM,Im«} plane. Black circles denote u\ roots, 
gray circles denote 114 roots. Vertical spacing between roots is i. 



In the thermodynamic limit, solutions of the above ABA are described by complexes 
of Bethe roots. Among those complexes are o„,A n ,ffi,®, which are the same complexes as 
in the mirror AdS5/CFT4 (see pSfl). In addition, the momentum-carrying roots u 4 and 
ug form two new types of complexes, which we call Odd-Even and Even-Odd. They were 
recently considered in [ 44 ] . Those complexes are real-centered strings of alternating u 4 and 
ui roots, adjacent roots being spaced by i. In the Odd-Even complex, the lowest root of 
the string on the complex plane is U4, while in the Even-Odd complex, the lowest root is 
u\ (see Fig. 2). The list of all complexes is given in the table below. 



On 


string of roots 


u 2 = u + ij, j = -^,...,v=Z 


\i 


pyramid 


u 3 = u + ij, j = -0=1,..., 2=1 
u 2 = u + ij, i = -2=2,..., 2=2 
ui=u + ij, j = -2=5,..., 2=3 


& 


single fermion root 


U\ = u 


8 


single fermion root 


u 3 = u 


•*n 


Odd-Even complex 


u 4 = u + ij when 6f n = l, j = -2=i, . . . , 2=1 
u\ = u + ij when 9f n = 1 


►n 


Even-Odd complex 


u 4 = u + ij when 9f n = 1, j = -2=1, . . . , 2=1 
ug = u + ij when 9f n = 1 



Here, «6l denotes the center of a complex, notation j = — 2k-=, . . . , 2— means that 



j takes the values —2=1, — 2=^, . . . , 2=^ 2=1^ an d 0' s we re defined in Q2.16| ). 

The energy (in our notation ip* A , with index A taking the values {o n ,e,®,A n , •«„,►„}) 
which corresponds to a complex is the sum of energies of the roots in a complex, and the 
same is true for momentum. We have p^ = p£ n = p* , = = e* , where 

KM - J ^ <(«) - f + MA) (-^-ttLtV (3-18) 



an 



while for other complexes, p* A and are zero. Note also that the only complexes with odd 
number of fermion roots are those denoted by ©, ®, <2n-i and *2n—l- Hence in our case the 



quantity hj^ in ( 3.11 ) has to be log(— 1) for these complexes and log(+l) otherwise. 
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Applying the fusion procedure to the mirror ABA equations 10 , we find that in our 
case, kernels Kab entering ( |3.11 ) are given by the table below. 



A\B 


°m 




8 


A m 


"*m 


*"m 








+K n -i 











© 


—Km-1 










R (oi) 


R (01) 


8 


—K m -i 








+K m -i 


77(01) 


77(01) 

/v lm 


A„ 









+i^n— l,m-l 


77(01) „(01) 
'v nm <-'n_2.m 


77(01) R (01) 
'v nm u n -2,m 


■*n 





R (10) 
D nl 


,v nl 


77(10) K (10) 

/Vnm i3 n,m-2 


_-rll 

/ nm 


-r x 

' nm 


►n 





R (10) 


7? (10) 
'Si 


77(10) R (10) 
/Vnm "n.m— 2 


' nm 


-r" 

' nm 



Some of those kernels are the same as in the AdS5/CFT4 case [24], and we list them in 
Appendix A. The new kernels are 

7-II = r(h) _i_ s — k\\ 

1 nm — LJ nm ' '-'nm ^nmi 
«( n )+ 5 -K 1 

'-'■nm t <->nm lentil! 



' nm 



(3.19) 
(3.20) 



where 



n-1 
2 



771—1 
2 



.Kf m (u,v) 
K nmi U ' V ) 



n — 1 r, m—1 

2 K_ 2 

n — 1 m — 1 

£ e ^-ni(i-^)(ct+cc). 



(3.21) 
(3.22) 



m—1 
2 



Let us also introduce the functions Ya, which will turn out to be the functions which 
enter the Y-system: 



Yo n ,Ycs, — , — — , — , — J> = {3^,34,3^)3^,3^,34^} 

*® *A„ I-*n I>„ 



(3.23) 



(We recall that 3^ 



PA 

pa 



•) 



We can write the TBA equations for the ground state in the following way: 

1 + l / Y O m , ,,,(01) 



log!® = +iT m _i * log 
log!® = -iT m _i * log 



1 + l/%„ 



+ ^ * log(l + y^)(l + rO + ivr (3.24) 



^ ) *io g (i + y^)(i + y. i 



log 3^ = -K n _i im _! * log(l + Y A J - #„_i®log(l + Y 9 ) 

+ (n^ + 2>m ) * io g (i + y«j(i + ro 

log y Qn = #n-l,m-l * log(l + l/lb m ) + tf n -l®]og(l + 3£) 

io g y^ = -zL P ; + rJL*iog(i + ^)+r r i*iog(i + ^ m ) + 
+ io g (i + + (7^ + b£2_ 2 ) * i°g(i + ^m) 

log = -»Lp* + rJL * iog(i + + r n i * i og (i + Y <m ) + mn 
+ <\ 0) ®iog(i + + + B£2_a) * io g (i + y Am ) 



27T 



(3.25) 
(3.26) 

(3.27) 
(3.28) 

(3.29) 



3 an important assumption here is monotonicity 
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where * denotes integration over the second variable, as in ( 3.11 ), Summation over the 
repeated index vn is assumed with ui J> 2 for is m and Om , and tti ^ 1 for < m , > m . 

Range of integration for fermions is limited to —2h < u < 2h. Notice that from 



( |3~24|) and ( |3~25| ) we can see that A is the analytical continuation of 1® across the cut u £ 



(— oo, — 2/i) U (2h, +oo). For the convolutions with fermions we introduce the convolutions 
<£> which should be understood in the sense of a B-cycle (see pl|), e.g. 

/2ft ^ _i_ y 

dvK n -i log — — fr , (3.30) 
-2/j i + i/ r© 

/•2ft. 

7^ n( %iog(i + y») = / d« n (n ^ io g (i + y») - £ (n0) io g (i + i/y©; 

./-2ft L 

Remarkably, the combination Bum + S nm , which is part of the kernels T^m an d Tnm, 
has only two branch cuts for each of the variables u and v. This follows from the integral 
representation 

S nm (u,v)+B^(u,v) = (3.31) 



oo „ 

£ / [b£ 0) (u, u; + ia/2) £ff («, - ia/2, v) + B^ (u, u, - ia/2) sff (w + fa/2, ,,) 

o=l ^ 



diu, 



which can be derived using the results obtained in [^4], ^] (see Appendix A). As a conse- 
quence, the functions Y^u) and Y+. n {u) (see ( 3.28[ ), ( |3.29|) ) should not have branch cuts 
for —in/2 < Imu < in/2. 

In the next section we will establish a relation between the above equations and the 
AdS 4 /CFT 3 Y-system. 

3.2 Y-system from TBA equations 

In this section we show that solutions of ground state TBA equations satisfy the AdS4/CFT3 
Y-system ( |1.2| )-(|i~4|) described in the introduction. This derivation of the Y-system is sim- 
ilar to the AdSs/CFT4 case [24]. First, we identify the Ya functions in TBA equations 



with the Y-system functions Y a<s . We set 

{Yo^Y^Y^J = {Yi,n,Y2,2,Yi,i^n,i}- (3-32) 

Let us introduce the discrete Laplacian operator 

AK n {u) = K n (u + i/2 - »0) + K n (u - i/2 + iO) - K n+1 (u) - K n ^(u). 



Following [E4J, we apply this operator to the l.h.s. of the TBA equations, acting on the 
free index n and the free variable u. The action of this Laplacian on some of our kernels 
has been computed in 

AK n (u) = 6 nA 5(u) (3.33) 

AK nm (v - u) = ATZ^(v,u) = 5 n>m+1 5(v - u) + <5 n)m _i<5(v - u) (3.34) 

ATl^(v, u) = Ak£$(v, u) = 5 n , m 5(v - u) (3.35) 

AB nm = 0, AS nm = . (3.36) 
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The new kernels Tnm and 7^ satisfy relations of a new type, which are not written in 
terms of the Laplacian: 



* I n.rn \ ^ ~~ l~~ 



i-iO 
2 ' 



-0 "^"^nmC^ 



i-iO 
2 ' 



(3.37) 



T^ m {u+^,v)+Tn m (,u- l -^f,v)-Tn+i, m (u, «)-7ii )m (u, «) = -<y n -i,in*(tt-«) • 

Using those identities, we obtain from the TBA equations a set of simpler equations for 
the functions Ya- This closely follows [24]. For example, applying the Laplacian to the 
l.h.s of (ET27D , we get 



log 



Yc^+iY^-i 



log(l + l/y „ +1 )(l + l/yo„_ 1 ) , n>2 



or, equivalently, 



logY+Y c r=log(l + ycw +1 )(l + ^o n _ 1 ) , n>2. 



(3.38) 
(3.39) 



For n = 2 we obtain 



log 



(i + y 8 )(i + y 



1 + I/Yb 



(3.40) 



Equations (3.24) and (3.26) can be treated in a similar way. We get an equation for Yg 

(i + YoXi + ^xi + hj 



log 

and also equations for Ys n : 



log- 



1 + V^A 2 



log 



YA n+ iiA rl _i 



(1 + yQCi + yQ 



(3.41) 



(3.42) 



log- 



Y A +Y: 



log 



(i + y ffi )(i + K t )(i + i^)y 



(i + Ya 3 )(i + y®) 
log y«y ffi + - * lo sd + y OU + Kv 



(3.43) 



Moreover, adding up Eqs. (|3.24| ), fl3.25|) we find that 

log Y®Y® = E^i- " fi im) * lo g(! + ^0(1 + 



(3.44) 



Therefore, in Eq. ( 3.43j ) all summands except the first one cancel, and that equation takes 
the compact form 

:i + y e )(i + y,)(i + y^) 



log ISPS 



log- 



(3.45) 



- - (i + i/y A3 )(i + i/y) ' 

Equations for Y^, Y+ n are obtained from ( |3.28 ), ( |3.29| ) in a similar way with the use of 
new identities ( 3.37] ), and they are precisely equations fll.2| )-(|L"4|) which were given in the 
introduction: 



1 + Ya„ 



i + y A „ 



iogy+y.: = io g . 



l + i/y^Xi + i/y^O 



, n > 1, 



(3.46) 
(3.47) 
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while for n = 1 



log Y^+Y^ =log— ^— — , log Y+Y^ =log- 



(1 + V^) 

3.3 Integral equations for excited states 



(1 + V^) 



(3.48) 



As we have shown above the equations ( 3.24| )-( |3.29| ) contain important structural infor- 
mation about the Y-system. However, those equations do not make much sense when 
understood literally since they describe the ground state which is protected by super- 
symmetry, and the Y-functions are degenerate in this case. There is a way to extend these 
equations to excited states, with the Y-functions becoming very nontrivial. For the case 
Y^ = Y*. a = Y# a , Uij = ugj similarly to [24] we propose, as a conjecture, the following set 
of equations: 



1"A V +/v',„- i * log 1 *yj° m + 2K { S * log(l + Y.J + 2 



loglffi 



-K, 



m—l 



* log 



1 + Ya 



log 



R 



(+) 



R 



(-) 



2^° m 1 )*log(l + YO-2 



log 



B 



(+) 



B 



(-) 



logY An = -K, 



n— l,m— 1 



* log(l + Y A , 



K„_i * log 



+ 2(^) + 4 _ 1 ) im )*iog(i + Y^; 



1 + 1® 
1 + l/Y m 



+ in (3.49) 

- m (3.50) 
(3.51) 



+ 2 



n-l 

2 



E ^ 



R { +\u + ik) 



+ 2 



n — 3 
2 



E ^ 



ra — 3 
2 



ff^V + zfc) 
sj _) (u + ^) 



1 + Y 8 



log Yw = K n - 1>m -i * log(l + 1/lcw) + #n-i * log , , 

J. + -L/iffi 

[— n] 

log Y^ = Jlog ^ - * log(l + 1/Ye) + <\ 0) * log(l + Y A 
+ (^)+^l 2 )*log(l + Y Am ) 



(3.52) 
(3.53) 



(2S nm - + B<#) * log(l + Y.J + 



2 



y~] log$ 4 (w + i/c) 



2 



+ i7rn 



where J = L -\- $4 is given by (|2~T7|) (with Bf,Bf in that expression replaced by 
unity) and the exact positions of the Bethe roots are determined by 



-1, j = l,...,K 4 . 



(3.54) 



The label "ph" here means that one should analytically continue the equation for Y^ to the 
physical sheet, like it was done for the first time in |29[| . The Bethe roots are additionally 
constrained by a condition imposed on total momentum (the trace cyclicity condition) . We 
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can write this constraint in a form similar to ( |1.5| ): 



00 f°° dn f)n mil (u) Ka 

a=l J -°° j=l j=l 

(3.55) 

(recall that the momentum p(u) was introduced in ( 3.13[) ), This expression can be simplified 
in our case, as Y^ = Y*. a = Y 9a and u^j = u%j. 

Note that in equations for excited states, in the terms without convolutions the branch 
x ph should be used for x(u4j) and x mir should be used for x(u) with u being the free 
variable. 

Strictly speaking these equations are only valid for some particular values of A and 
configurations of roots. In other cases the equations may require some modification. This 
question is usually subjected to case-by-case study (see e.g. []55], 54, 53]). 

In general the procedure is the following - one can start from a sufficiently large L or 
small A where the terms with log(l + Y <n ) are irrelevant and the asymptotic solution of [19] 
should be a good approximation. The condition fl3.54|) can be discarded for a while, and 
one should find such a configuration of the roots U4 (usually they are sufficiently close 
to the origin in this case) that the asymptotic solution satisfies the equations for excited 
states we proposed above. After that the equations should be analytically continued in 
L, A and u^j. 

This procedure in general is rather complicated, however our experience with the 
Konishi operator in AdSs/CFT4 [29] tells us that one can probably use the equations 
above as they are from A = to very large A's. At the same time the Y-system functional 
equations are not affected by these modifications and they are more suitable for the strong 
coupling analysis [|32]]. Moreover, they are not restricted to the "s/(2)" subsector. 

The possibility that some singularities could collide with the integration contours and 
modify the equations when some parameters (such as the coupling) are changed was studied 



in detail in |^, 54, 53]. For AdS/CFT, this issue was mentioned in |29j| , and following that 
proposal, such a possibility was explored in [|58| for AdSs/CFT^ 11 



4. Solution of the Y-system in the scaling limit 

In this section we obtain a solution of the AdS4/CFT3 Y-system in the strong coupling 
scaling limit, considering the sl{2) subsector. In this case, the Y-functions which correspond 
to the momentum-carrying roots are equal. We show that the spectrum obtained from the 
Y-system is in complete agreement with the results from quasiclassical string theory. 

11 In |^8| an attempt was also made to estimate the "critical" values of the 't Hooft coupling - values 
for which the equations for excited states should be modified by extra terms. The result from [^SJ is 
A cr iticai — 774. The method used in that work is based on the asymptotic solution of the Y-system. 
The asymptotic solution works perfectly for very small and very large values of the coupling, however it very 
badly approximates the exact Y-functions for A ~ 700. Thus the only reasonable estimate at the moment 
for the critical value is A > 700, from the results of pflft , where no singularity was found in numerical studies 
of the TBA equations in the range < A < 700. 
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4.1 Y-system equations in the scaling limit 

In the scaling limit the Y-system simplifies in several important ways. In this section and 
section 5 we use rescaled rapidities z = 2 h(\) (similarly to @|), and since h(X) — > oo, we 
can neglect shifts in the arguments in the l.h.s. of the Y-system equations 12 . Hence with 
1/h 2 precision the Y-system becomes a set of algebraic, instead of functional, equations. 
Moreover, for the si (2) subsector Y^ = Y> a . Also, only z^ k and zg k Bethe roots are 
introduced (see (|l.6|) ), and they coincide pairwise: z^k = z^ k . Denoting Y 9a = Y^ we get 
three infinite series of equations 

*3, = + + y q ,_ I ) , « = 3,4 (4.1) 

^' (HVUM/^) '-"- (4 - 2) 

Y 2 = ; , r , a = 2, 3, ... , (4.3) 

* (l + l/Y. a+1 )(l + l/Y. a _ 1 ) { } 



plus four more equations 



Y 2 (1 + Ye)(l + YQ 2 

(i + i/Ya 3 )(i + i/y ) ' { -> 

2 _ (1 + Ya 3 )(l + Y 8 ) 
Y °> ~ (1 + 1/Y ffi ) ' (45) 

2 (1 + YcXi + y.j 2 



2 = (i + y®, 
* (l + i/V 



Together with the Y-system, we have to solve the non-local equation 

log Y®Y® = 2 ^ (7^ } - fig?) * log(l + Y.J + 2 log * +) , (4.8) 

m=l -^4 -^4 

which can be obtained by adding up ( 3.5CQ and ( |3.49|) (it corresponds to the gray node in 
Fig. 1, right). Introducing the following notation 

1 Mfe 1 x 2 

Gk{x) = -Y J — ^tt, ^ = 4,4 (4.9) 

n x x k - i 

f k (z) = exp (- »G fc (*(*))) , /*(*) = exp ( + iG k {l/x{z))) , (4.10) 

where the mirror branch of x is used for x{z) and the physical branch for x k j (this choice 
of branches is used by default in sections 4 and 5), following [32] we can write the non-local 
equation in the form 

1 oo 



n=l 



12 This simplification of the Y-system involves certain subtleties, as the shifts in the argument of the 
Y-functions cannot be neglected close to the branch cuts. This issue can be treated in our case in the same 
way as in (B2[. 
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where 

F = Y«&, f(z) = fi(z), f(z) = f!(z). (4.12) 
Note that, similarly to |52], the Bethe roots have to satisfy the constraint 



where 



Q 1 = 2irm + 0(l/h), m G Z , (4.13) 



n, v - ^ • >: - io, ^ - y — p— • >: ,, y ,, • (4.14) 



3=1 " i=1 * -4,j j=l /l ( X 4j !) j=1 ^4,j 

This condition reflects the cyclicity symmetry of single trace operators. Its consistency with 
the other equations remains to be checked, and we assume that equation to be satisfied. 

4.2 Asymptotics of Y-functions 

The Y-system equations should be supplemented by boundary conditions on the functions 
Y ajS , i.e. by their large a,s asymptotics. In our case, similarly to AdS5/CFT 4 (see |B3]), 
we demand that Y/\ a and Ym a have the same asymptotics as the L — > oo solution, which was 
constructed in section 2. As for the functions Yd s , we demand that their large s asymptotics 
is polynomial in s, which is true for the L — > oo solution as well. 

It is straightforward to show that the expressions for Y-functions from section 2 can 
be recast in the following form: 

f(f- 1)2 fa _ fff_ ^2 fa 

Y 0s (.) = (s - Af - 1 , Y Aa = {ST : +1 _^ {STa . 1 _ 1) , (4-16) 



where 



1 f f(f-l) 2 f 



while /, / are given by ( 4.12 ). As for real z the quantity f(z)/f(z) is a pure phase, to 
investigate the large a, s limit we consider Y-functions of shifted argument z — iO. We have 
then \f(z — iO)/ f(z — iO) > 1 and we get 

nm bgY^-zO) =log / 

a— >oo a J 

Similarly, 

lim logY * ( *- i0) =log(-A/) ■ (4.20) 

a— >oo Q 

Conditions ( 4.19|) , ( [1.20 ) are the boundary conditions which we impose on the functions 
Y\ a ,Y 9a for finite L at strong coupling. 
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4.3 Solution in upper and right wings 

In this section, we solve the Y-system partially, expressing the upper wing functions Y/\ a , 
Y% a and the right wing functions Yq 3 in terms of only three yet unknown functions. 

Using the analogy between our Y-system and the one considered in [32], it can be 
shown that the functions Y\ a ,Y 9a can be constructed in the following way: 

l + Y Aa = ^ , 1 + Y. a = ^ , (4.21) 

-^0+1,0-^0-1,0 -^0+1,1-^0-1,1 

where the set of functions T a>s , which is the general solution of the Hirota equation in the 
vertical strip, was found in [32]. Those functions are: 

T a , 2 = 1 (4.22) 

rp Vm ( Vl fg a 1 \ V2 ( q g i 

" (yi-toKvm-i) \ m + s^)~^F~i V 2V2 + 

Ta,o = (Tlx ~ T a+1>1 T a _ ltl )/T aj2 . (4.24) 



with ux, y2, Si and S2 being arbitrary parameters. The functions li a ,lA a , given by ( 4.21| ), 
satisfy the Y-system equations (|4.2| ) and ( fO| ) for arbitrary yi(z), 2/2(2), S±(z) and S2(z). 
The asymptotic conditions ( 4.1S| ),( fOo| ) fix y\ and y2'- 

Vi = ~T > V* = ■ (4.25) 

The general solution of ( f4.1| ) with polynomial large s asymptotics is (see [32]) 

Y 0s (z) = (s - A(z)f -1 , (4.26) 

with arbitrary A(z). Thus, the solution of our Y-system in the upper and right wings is 
expressed in terms of three unknown functions Si(z), S 2 (z) and A(z). 

4.4 Matching wings 

By now, we have constructed Yo n ,Y\ n and Y^ for all n in terms of A(z), S\(z) and S 2 (z). 
To find those three functions, as well as Y® and i®, we have to solve the five remaining 
equations (gg),(gg),(gjg),(gjg),( glg) . Excluding Y m and Y®, we get: 



Y^(1 + 1/Y A3 ) _ F 
(1 + Y^y (A - l) 2 

1 + VjAg 1 ( (A -I)* ^ 2 
(1 + Y % ) 2 A 2 \ F 



(4.27) 
(4.28) 



^(l + l/^) = - ( ^_ 1 ( ) }! 1)2 1) , (4-29) 
1 °° 



n=l 
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The r.h.s. of the four equations above depends only on F{z) and A(z), and they can 
be solved perturbatively in A, like analogous equations in j32|. Namely, we find several 
terms in the expansion of unknown functions in powers of A, notice a simple relation 
between consecutive terms, and sum up the series in A assuming this relation to hold for 
all terms 13 . It is then easy to check that the functions obtained in this way are indeed 
solutions of ( pL27p , flOf ), ( ^29| ). The result is: 



(1 + A)((l - / 4 / 4 )(l - A 2 / 4 / 4 ) ~ AC/4 - h?)((l + / 4 / 4 )(l + A 2 / 4 / 4 ) - A(/ 4 + hf) 
(-1 + A)(/ - 1)(/ - 1)(AV - 1)(A 2 / - 1) 

(4.31) 

(/-!)/( A 2 / -1) 

Sl ~ /(A 2 / -!)(/--!) (432) 
_ (/-i)(/-i) 

62 " //-(A 2 /-1)(A 2 /-1) (433) 

„ (A/-1) 2 (A/- l) 2 , . 

F = — ; \/ - ; (4.34 
(A -1)4// 1 ) 



Putting those functions into the expressions for the Y- functions ( 4.21| ), ( |4.26j ), we obtain 



all the Y a>s in terms of /, / and A (using (|4.7|) to find Y® and then getting Y$ from F = Y®Y®). 



Recalling the definitions fl£l!j ), (pl8|) of /, / and A, we see that we have found all 



the Y- functions in terms of the Bethe roots. We present our solution of the Y-system in 
Mathematica form in Appendix C. 

4.5 The spectrum from Y-system 



Here we repeat the arguments of ]32j to find the equation for the displacement of Bethe roots 
due to the finite size effects at strong coupling. In this section we assume Y* n = Y^ = Y^ ii 
(the situation where this is not the case is considered in the next section). We again start 
from the TBA equation for the momentum-carrying node 

logY % =T lm *log(l + Y^J + n^^log(l + Y^ + n^m m -i*log(l + Y Am ) + i^ , (4.35) 

where 7i m = 25 nm — IZnm + Bnm and $ represents extra potentials in the TBA equations 



for the excited states. The only difference with |32|] is absence of 2's in front of the second 
and third terms. Thus we can use the same trick as in [^] to get the expression for Y^ in 
physical kinematics: 

log |£ = 7^ * log(l + *J + 7^°)p^'© log ( j±_|) (4.36) 
+ K^™m m -i * log Q^|jr) + K ^z k - i- h ) * log (y^) • 



13 For A = 0, there are several solutions of Eqs. (4.27)-(4.29). We choose the one consistent with the 
asymptotic solution of Y-system. 
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Now we simply have to expand the kernels at large h and substitute Y's. Let us denote 

x 2 d 7 1 , x d 7 1 



r(x, z) 



, u(x,z) 



x 2 — 1 2irh x — x(z) ' v ' ; x i — i 2nh x 2 (z) — 1 
We rearrange the terms in ( [4,36 ) to evaluate the following "magic" products 14 

m=l 



D +M- 



1+jg TT 

AA U + W 



1 + ^° 

8 m=2 

1 + 1/Y£ t r / l + V 



1 + 1/% 



n i 



m=2 



(A - 1)6(A + 1)2(A/ + 1)2(A/ + 1) 2 (A2// - I) 2 

m fr 1 (/A + l) 2 (//A 2 -l) 

AA (l + ^) 2m (/A 2 -l) 2 

m— 2 / Ja2 i\2 



(/A 2 - 1)> 



(/A + l) 2 (//A 2 -l) 



(4.37) 



We get the following corrected Bethe equation for the 5I2 sector 



L M 



a n 



41 



J= i x k 



—a {z k ,Zj) 



x exp 



^ 1 - 1 /( x k x j) 
(r(x kl z)M + - r(l/x k , z)M- + u(x k , z)M Q S jdz 



and the equation for the energy at strong coupling is 

-1 dz 
-1 



E 



M 2,i 



i=l 



X? 



1 



(4.38) 



(4.39) 



The extra factor of 2 in the denominator under the integral is due to the single magnon 
dispersion relation, which includes an extra 1/2 compared to AdSs/CFT^ 

4.6 Non-symmetric strong coupling solution 

In this section we present a simple strong coupling solution with Y+^Y*.. It can be used 
as a test of the new structure of the Y-system which we proposed in the introduction. 

We consider the limit where the massive nodes are completely decoupled from the rest 
of the system. We solve the following infinite set of equations 

1 



;i + i/y^ +1 )(i + i/i; n „ 1 ; 
1 



n > 1 
n > 1. 



(i + i/Y. n+1 )(i + i/i^_ 1 ; 

The explicit general solution of this system with two parameters a and (3 is 

(a 2 -l)V 



(4.40) 
(4.41) 



Y, 



(a 2 -l)V 
(o n - 1) (a+ (a n + 2 - 1) 



n is odd 



n is even 



(4.42) 



14 To compute these products we again use the z — iO prescription to ensure their convergence. This 
prescription is inherited from the TBA equation for excited states where the integration should go slightly 
below the real axis. 
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Figure 3: Elementary excitations of the string in AdS4 x 



and Y+ n is obtained from by replacing /3 — >■ 1//3. We can easily compute .Mo for this 
solution, 

c M _ (aP + l)(a + P) 
(a 2 -l) 2 /? 

and by matching with the asymptotic solution we identify 



(4.43) 



so that we get 



Ai = a , A 2 = A 3 = A 4 = , A 5 = -(3 . 
„Mo_ (A 5 -A 1 )(1-A 1 A 5 ) 



(A?-1) 2 A 5 



(4.44) 
(4.45) 



5. One-loop strong coupling quasi-classical string spectrum 



In this section we briefly describe the construction of [32], applied for the ABJM model. 
The algebraic curve described in can be used to compute the one-loop correction for a 
generic finite gap classical string state by computing the spectrum of fluctuations around 
a given solution. We assume that the one-loop shift computed from the algebraic curve 
agrees with the strong coupling expansion of ABA in the limit L/h S> 1 . This assumption 



was explicitly verified for the folded string in [59]. The general proof like in [60] is still 
missing. 

There is yet another way to compute the one-loop shift directly from the world-sheet 
action which is similar to the algebraic curve computation. Whereas for the folded string 
both computations give the same excitation frequencies 15 [61, 59], for the circular string 
a negative result was obtained in [62]. Recently it was shown that one should be more 



15 similar analysys for the giant magnon was done in |3q, 
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carefull with the periodicity of the fermionic fields in the world-sheet approach and the 
corrected derivation leads to agreement with the algebraic curve frequencies Q so we 
assume all approaches to be consistent with each other. 

The pattern of excitations in the ABJM theory is quite different from that of AdSs x S 5 . 
The string in AdS4 x CP 3 has 8 bosonic (3 modes of AdS4 and 5 of CP 3 ) and 8 fermionic 
excitations. They are divided into heavy and light modes (see Fig.|3|). The dispersion 
relations for the heavy and light modes differ by a factor of two. We will see that this 
complicated structure of heavy and light fluctuations is captured by the Y-system. As 
usual the one- loop shift is given by a sum over the fluctuation energies |64| . In the algebraic 
curve language the fluctuations are the small cuts (i.e. poles) connecting different sheets 
of the algebraic curve The poles could be placed only in certain special positions x % £ 
given by 

«(«n ) - Qj(?%) = 27m • 
The quasiclassical Bohr-Sommerfeld quantization condition constrains the minimal residue 
of the pole. Insertion of the pole results in displacement of the other singularities. Moreover 
the pole by itself carries the energy co(xn) for the light mode and 2tu(xn) for the heavy 
mode where 

x z — 1 



Following [32] we first compute this second part of the one-loop shift which does not take 



into account the back-reaction of the fluctuation on the large cuts. Then we have 

sei 1oop = \ e Ec-^m^') + E E (-i^otf) (5-3) 

n light n heavy 

where Fij = 1 for bosonic modes and —1 for fermionic. The modes are 

(ij) = (4,5), (4,6), (3,5), (3,6) , light bosonic modes (5.4) 

(ij) = (3, 7), (2, 9), (1,9), (1,10) , heavy bosonic modes (5.5) 

= (2, 5), (2, 6), (1, 5), (1, 6) , light fermionic modes (5-6) 

= (2,7), (1,7), (2,8), (1,8) , heavy fermionic modes (5.7) 

Rewriting the sum ( |5.3[) as an integral we get 

SEtioop = E f ^»(*)Wo (5-8) 

where the integration goes over the upper half of the unit circle \x\ = 1 and we denote 

e M) = JJ(1 — e ~ ipi+ip ^) Fi ^> JJ (1 — e ~ iPi+ip ') 2Fi ' (5.9) 

light heavy 

Note that d x Afo is constructed to have the residue ±1 (±2) exactly at the position of the 
light (heavy) mode Xn ■ More explicitly we can write 

A/- (Aj - 1) 2 (A| - 1) 2 (A!A 2 - 1) 2 (A 3 A 4 - 1) 2 (A 3 - A 5 )(A 4 - A 5 )(A 3 A 5 - 1)(A 4 A 5 - 1) 

" (AiA 3 - 1) 2 (A 2 A 3 - l) 2 (AiA 4 - 1) 2 (A 2 A 4 - 1) 2 (A! - A 5 )(A 2 - AsXA^s - 1)(A 2 A 5 - 1) 

(5.10) 
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where 

\ a = e~ iqa . (5.11) 
For s/2 sector |12[ | there are only cuts connecting 2nd and 9th sheets so that 

A 5 = 1 , A3 = A4 = A , A 2 = A/ , A 2 = Af (5.12) 
and ( |5.10| ) simplifies to 

Nn _ (A - 1) 6 (A + 1) 2 (A/ + 1) 2 (A/ + 1) 2 (A 2 //- l) 2 ir. 1 o\ 

(A 2 /-1) 4 (A 2 /-1) 4 ' { ' 

where we recognize Mo = — Mio and ( |5.8| ) coincides precisely with the second term in ( 4.39| )! 

Then one should also take into account the back-reaction of the fluctuations. As it is 
explained in detail in |3^] for that one should work with the modified Bethe equations. We 
need to consider only the fluctuations touching those sheets where the macroscopic cuts 
are located. Hence, one of those sheets has to be the 2nd or the 9th sheet. Computing the 



r.h.s of (|5.9|) with that restriction imposed, we get, similarly to [32] 



e M + = (AiA 2 - 1)(A 2 - l) 2 

(A2A3 - 1)(A 2 A 4 - 1)(A 2 - A 5 )(A 2 - I/A5) 1 • ' 

e N. = (A 2 - 1) 2 (A!A 2 - 1) 

(A1A3 - l)(AiA 4 - l)(Ai - A 5 )(Ai - I/A5) 1 • ' 



and from ( f>. 12 ) we get 



x+ _ (A/ + 1) 2 (A 2 //-1) „_ (A/+1) 2 (A 2 //-1) 

(A2/-1) 2 ' (A 2 /-l) 2 [ } 

and we recognize exactly the same structures ( 4.37] ) we got from solving the Y-system! 



Finally by putting A2 = A3 = A4 = we obtain 



(A 1 -A 5 )(A 1 A 5 -1) { n 

which is again precisely the quantity e _-A/! ° obtained for this sector in ( 4.45| )! 



We see that the nontrivial pattern of the fluctuations is reflected in the Y-system thus 
providing a direct link with the worldsheet theory. This is also a deep test of the structure 
of the Y-system equations we proposed. 

6. Conclusion 



In this paper we refined the Y-system for the ABJM theory which was conjectured in [19]. 
We derived it directly through the TBA approach and then made several highly nontrivial 
tests at strong coupling. In particular we constructed the general 5I2 solution for the new 
Y-system in the scaling limit, and also made a test for a subsector where the difference 
between the old and the new Y-systems is crucial. 

We also constructed the general asymptotic solution of the Y-system for arbitrary 
excited states. It can be used, in particular, for the weak coupling tests of the conjecture 
and as an initial configuration for numerical iterative solutions. 



Note added While we were working on the strong coupling solution, the paper [65] 
appeared, with a similar Y-system and vacuum TBA equations. 
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A. Notation and kernels 



We use the following notation for kernels in TBA equations: 

K n (u,v) 



where 



and 



1 d u — v + in/2 
- — In ■ 



27ri dv u — v — in/2 

Tn— 1 n — 1 

2 2 

K n>m (u,v) = ^ Yj K 2j+2k+2{u,v), 
Snm(u,v) = J-^-loga BE s(x^(u),X^ n \u),X^(v),X^(v)) 



2iri dv 



S nm (u,v) = S nm (u,v) + ^-V {m) (v) 



n—1 m—1 
2 2 



1 d j b(u + ia/2 + ij,v — ib/2 + ik) 
^ ± 2iri dv ° g b(u - ia/2 + ij, v + ib/2 + ik) 



n—1 i m 



n—1 m—1 

2 2 



1 d j r{u + ia/2 + ij,v — ib/2 + ik) 
2—* ^ i 27ri dv °^ r(u — ia/2 + ij, v + ib/2 + ik) 



■ _ n—1 h.__m- 

3- K ~ — 



,. . l/x mir (u) -x mh (v) . . x miT (u)-x m[T (v) 
b{u,v) = . , r(u,v) = - 



^x miT (v) 



x mh (v + ia/2) 



w 2?r & x mir (v - ia/2) 



(A.l) 

(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 

(A.7) 
(A.8) 



To derive fljDfij ) we use the following integral representation |24|, |25f : 
25 nm (u, u) - (u, v) + (u, v) = -K njWl (u - v) 

oo „ 

-2 £ / [«£° («, to + fa/2) Bff (to - ia/2, «) + fig" (ti, w - ia/2) sff (to + ia/2, «) 



a=l 



(A.9) 



Equation ( p. 31 ) follows due to the identity 1Znm (u, v) + Sim (u, v) = K njrn (u — v) 
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B. Fermionic duality transformation and su(2) 

We can transform a set of Bethe equations into an equivalent one by application of the 
fermionic duality. This follows || closely. We construct the polynomial 



K 4 


K- 4 


K 2 












3=1 


3=1 


3=1 




K A 












~ X ^')II( X ~ X ^3^ X 




3=1 


3=1 


3=1 





(B.l) 



from the Bethe equations of |12| (given in section 2.1 of the present work) for the fermionic 
roots u% and u 3 . We see that this polynomial has zeros for x = l/x\j and x = x 3 j. 
Denoting the remaining zeros of this polynomial by x, we get 



K 3 Ki Ka Ki 

P(X) = C*[[(x- X 3 j) 1[(X- 1/X hj ) 1[(X- X 3J ) ]J(X - l/ Xld ) 
3=1 3=1 3=1 3=1 

(where C is a constant) or in our notation (with R = B = B4B4) 

K 4 



(B.2) 



P(x) = CR^R^Bi = R { - ] Qt - R (+) Q 2 



x\ 
h) 



K 2 



n 

3=1 



R^B^RtBt 



and 



then 



P(x + 

PjxF) R 3 B^R-B^ 
P(l/x-) _ B 3 R^B^R7 



P(l/x+) B+RfB±R± 

+ \K 2 R (+)+ B -b-B+B. 



K2 fi(-)+Q++ - i?(+)+Q 2 
R(-)-Q 2 - RM-Q2- 

I<2 B^-)-Q 2 - BM-Q 2 - 
B(-)+Q++ - B(+)+Q 2 



a-l 
2 



R(-)- B+ B+ B% B 3 



U u ) = II /(« + 



in) 



T a ,l(u\{ui t j},{U3,j}) = fa(u)T ha (u\{u ltj },{u 3 j}) 



(B.3) 
(B.4) 

(B.5) 

(B.6) 

(B.7) 
(B.8) 



T hs (u\{ui d },{u 3J }) = f s (u)T Stl (u\{u lt j},{u 3;j }) 
Here, the bar means complex conjugation in the physical sense, i.e. the replacement 

#)±^^)T } #)±^bWt. (b.9) 
Notice that x is not inverting under this conjugation. 

/ , , \ L-K 2 



Y^ 



T 0) i [| $4 (ii + in)<lj na (M + m), 



(B.10) 



L-K 2 



a-l 
2 



T a>1 ]~| $4" a (u + in)^" a (u + m) (B.H) 
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As in section 2, the factors <&4(u) and are constructed in such a way that the ABA 

equations for the momentum carrying nodes are given by Y^ h (u4 j) = — 1 and Y^^uj ■) = 
— 1. We have 



r (+)+d(-)-5-5+ B (+>+ R ( ~>~ B~B + 

(B.12) 



where Bi,Ri are defined similarly to (|2.8|), (|2.9|), with xjj replaced by xij. 



C. Explicit expressions for Y- functions 

Here we present the solution of the Y-system in the scaling limit for the si (2) sector. This 
solution was obtained in section 4, and below we give its explicit form, which can be used 
in the Mathematica system. We denote d = A, g = / 4 , gb = f^, Ym[a_] = Y 9a , Yp[a_] = 
Y& a , Yb[s_] = Yq 3 , and the Y- functions are given by the following code: 

sb={ 

A-> C(l+d)(l-d g"2-g gb+2 d g gb-d'2 g gb-d gb~2+d"2 g~2 gb~2) 
(1-d g~2+g gb-2 d g gb+d~2 g gb-d gb"2+d-2 g~2 gb~2))/ 
((-1+d) (-1+g) (1+g) (-1+d g) (1+d g) (-1+gb) (1+gb) (-1+d gb) (1+d gb)) , 

51- >((-l+g)(l+g)gb~2(-l+d gb)(l+d gb))/(g~2(-l+d g)(l+d g) (-1+gb) (1+gb)), 

52- >((-l+g) (1+g) (-1+gb) (l+gb))/(g~2(-l+d g)(l+d g)gb"2(-l+d gb)(l+d gb)), 
P-> ((-1+d g~2)~2(-l+d gb-2)-2)/((-l+d)-4g-2gb-2) , 

T2->d g gb, Tl->-(g/gb)}; 

Ym[a_]=-1+(S2 T2~ (1+a) (-1+T2~2) -Sl~2 S2 Tl~ (4+2a)T2~ (1+a) (-1+T2~2)+ 

SI Tl~(l+a) (-1+T1-2) (-1+S2-2 T2~ (4+2a) ) ) ~2/ 

((-S2 T2~a (-1+T2~2)+S1~2 S2 Tl~(2+2 a)T2"a (-1+T2"2)- 

Sl Tl~a (-1+T1"2) (-1+S2"2 T2~(2+2 a)))(-S2 T2~ (2+a) (-1+T2"2) + 

Sl~2 S2 Tl"(6+2 a) T2~ (2+a) (-1+T2~2) - 

SI Tl"(2+a) (-l+Tl'2) (-1+S2-2 T2~(6+2 a))))/.sb; 



Yp[a_]=((Sl Tl~(4+a) T2)/(-l+Tl~2)+(Tl~-a T2)/(S1-S1 Tl~2)- 
(Tl(T2~-a-S2~2 T2~ (4+a) ) ) / (S2 - S2 T2"2) ) "2/ (- ( (SI Tl"(4+a)T2)/ 
(-l+Tl~2)+(Tl~-a T2)/(S1 - SI Tl~2) - (Tl (T2~-a-S2~2 T2"(4+a)))/ 
(S2-S2 T2"2))-2+(Tl-(-2 a)T2"(-2 a) (T2"2-S2-2 T2"(4+2 a)- 
2 SI S2 Tl"(2+a)T2"(2+a) (-l+T2"2)+2 SI S2 Tl" (4+a) T2~ (2+a) (-1+T2~2) + 
Sl~2 Tl-(6+2 a)T2"2(-l+S2~2 T2"(2+2 a))+2 Tl T2 (-1+S2"2 T2"(4+2 a))- 
2 Sl~2 Tl"(5+2 a)T2(-l+S2~2 T2"(4+2 a))+Tl~2 (1-S2~2 T2~(6+2 a))+ 
Sl~2 Tl"(4+2 a)(-l+S2~2 T2"(6+2 a) ) ) ~2) / (Sl~2 S2"2 (-1+T1~2)~2 
(-1+T2~2) "2 (T2+T1~2 T2-T1 (1+T2"2) ) "2) ) / . sb; 



Yb[s_]=(s-A)-2-l/.sb; 
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Y11=(-1+(S1 S2(-1+T1~2) (-1+T2~2) (T2~2-S2~2 T2~6-2 SI S2 Tl~3 T2"3 (-1+T2~2)+ 

2 SI S2 Tl~5 T2~3(-1+T2~2)+S1~2 Tl~8 T2~2(-l+S2~2 T2~4)+ 

2 Tl T2(-l+S2~2 T2~6)-2 Sl~2 Tl~7 T2(-l+S2~2 T2~6)+T1~2 (1-S2~2 T2"8)+ 

Sl~2 Tl~6 (-1+S2-2 T2~8))~2)/((S1~2 S2 Tl~4 T2 (-1+T2~2)+S2 (T2-T2~3)- 

Sl Tl (-1+T1-2) (-1+S2-2 T2~4))~2 (T2~2-S2~2 T2~8-2 SI S2 Tl~4 T2~4 (-1+T2~2)+ 

2 SI S2 Tl~6 T2~4 (-1+T2~2)+S1~2 Tl~10 T2~2 (-1+S2~2 T2~6)+2 Tl T2 (-1+S2~2 T2~8)- 

2 Sl~2 Tl~9 T2 (-1+S2"2 T2~8)+T1~2 (1-S2~2 T2~10)+S1~2 Tl~8 (-1+S2~2 T2~ 10) ) ) ) / . sb ; 

Y22=(P/Yll)/.sb; 
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